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Abstract. We have combined infrared data with H/, H2 and H// surveys in order to spatially decompose the 
observed dust emission into components associated with different phases of the gas. An inversion technique is 
applied. For the decomposition, we use the IRAS 60 and 100 /im bands, the DIRBE 140 and 240 fim bands, as 
well as Archeops 850 and 2096 fim wavelengths. In addition, we apply the decomposition to all five WMAP bands. 
We obtain longitude and latitude profiles for each wavelength and for each gas component in carefully selected 
Galactic radius bins. We also derive emissivity coefficients for dust in atomic, molecular and ionized gas in each 
of the bins. The H/ emissivity appears to decrease with increasing Galactic radius indicating that dust associated 
with atomic gas is heated by the ambient interstellar radiation field (ISRF). By contrast, we find evidence that 
dust mixed with molecular clouds is significantly heated by O/B stars still embedded in their progenitor clouds. 
By assuming a modified black-body with emissivity law A"'^'^, we also derive the radial distribution of temperature 
for each phase of the gas. All of the WMAP bands except W appear to be dominated by emission from something 
other than normal dust, most likely a mixture of thermal bremstrahlung from diffuse ionized gas, synchrotron 
emission and spinning dust. Furthermore, we find indications of an emissivity excess at long wavelengths (A > 
850 fim) in the outer Galaxy (R > 8.9 kpc). This suggests either the existence of a very cold dust component in 
the outer Galaxy or a temperature dependence of the spectral emissivity index. Finally, it is shown that ~ 80% of 
the total FIR luminosity is produced by dust associated with atomic hydrogen, in agreement with earlier findings 
by Sodroski et al. (1997). 

The work presented here has been carried out as part of the development of analysis tools for the planned European 
Space Agency (ESA) Planck mission. 
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1. Introduction 

The release (Neugebauer et al. 1984) of the Infrared 
Astronomical Satellite (IRAS) maps (at 12, 25, 60 and 
100 /im) provided a new perspective on the Galactic 
interstellar medium (ISM) by revealing diffuse infrared 
emission over almost the entire sky. This finding was con- 
firmed (Hauser 1993) by the Diffuse Infrared Background 
Explorer (DIRBE) experiment (1.25 to 240 /im) on board 
the COBE satellite. Even in the lowest 12-/im IRAS 
band, stellar radiation contributes only a small fraction 
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(~ 8%) of the emission (Boulanger & Perault 1988), 
the rest is from ISM. This fact is also corroborated by 
the tight correlations found by Boulanger et al. 1996 at 
high Galactic latitudes, i.e. |&| > 20°, between the 100 
fim emission and the 21-cm emission, the latter being a 
faithful tracer of HI. By using a fixed spectral emissivity 
index /3=2, these authors derive a temperature for dust 
associated with HI of 17 K. A similar correlation has 
been found for high-latitude molecular clouds, their CO 
emission being tightly correlated with the 100 fim inten- 
sity distribution (e.g. Boulanger et al. 1996). Likewise, 
Lagache et al. (2000) detect dust emission from the Warm 
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Ionized Medium for b < -30°, b > 25°, with physical 
properties (temperature and emissivity) similar to dust 
associated with atomic hydrogen. Remarkably, most of 
this dust emission is not contributed by discrete H/j 
regions but by the diffuse gas. 

As one can notice, all the results mentioned above 
have been obtained at high latitudes. In fact, the Galactic 
plane in in general much more complicated to investigate. 
This is because at high latitudes each line of sight traces 
a relatively short path through the Galaxy and this 
makes it possible to probe each distinct phase of the 
gas separately. On the contrary, in the Galactic plane 
the lines of sight are much longer and the phases of 
the gas heavily blended. A solution to this problem 
has been provided by the maximum- likely hood method 
worked out by Bloemen et al. (1986). This method allows 
one to spatially decompose the total infrared emission 
into components associated with atomic, molecular and 
ionized gas in different bins of Galactic radiu. A similar 
technique is used in this paper and will be described in 
Section 2. 

Another important observational fact derived from 
IRAS data is that most of the IR luminosity falls at A 
> 60 /xm. In particular, ~ 50% of the IRAS luminosity 
is observed at A > 100 /im and only one sixth in each 
of the other bands. According to a variety of dust grain 
models (e.g. Mathis et al. 1977, Draine & Lee 1984, 
Draine & Anderson 1985; Weiland et al. 1986; Desert et 
al. 1990; Li & Greenberg 1997; Dwek et al. 1997; Draine 
& Li 2001; Zubko, Dwek & Arendt 2004) the 12 fim 
emission is dominated by PAHs (Polycyclic Aromatic 
Hydrocarbons) which are characterized by a size of order 
of a few nanometers, while at 100 fiia the emission is 
mainly due to large grains, typically in the size range 
10-20 nm to 0.1 /xm. Intermediate wavelengths (25-60 
/im) are instead the domain of VSGs (Very Small Grains), 
in the nanometer range. Therefore, IRAS data indicate 
that most of the emission is imdisputedly due to large 
grains. IRAS has also revealed a wide range of infrared 
dust colors from 12 to 100 /im (Boulanger et al. 1990). 
This can be explained through variations in the dust 
size distribution. However, the same phenomenon has 
also been observed at submillimeter wavelengths (200, 
260, 360 and 580 iim) by the PRONAOS balloon-borne 
experiment (Bernard et al. 1999, Stepnik et al. 2001) and 
it has been interpreted as due to changes in the physi- 
cal properties of the large grains themselves, i.e. grain 
composition and temperature as well as optical properties. 

The present work is motivated by several factors. In 
recent years there has been an unprecedented number 
of new centimeter and sub-millimeter experiments pro- 
ducing a large amount of data. Among these are the 
Archeops (Bcnoit et al. 2002) ballon-borne experiment 
and the WMAP ( Wilkinson Microwave Anisotropy Probe) 
satellite (Bennett et al. 2003a). Such experiments have 



provided maps of the sky in the frequency range 23 to 
545 GHz (i.e. ~ 10 mm down to 550 /xm) at an angular 
resolution, depending on frequency, from ^ 50 down 
to ~ 10 arcmin. Since most of the infrared emission is 
emitted at A > 60-100 /xm, the combination of these maps 
with existing IRAS and DIRBE maps now allows one 
to trace completely the dust spectrum in the frequency 
range where most of the emission is concentrated. At 
the same time, new all-sky data have also been released 
for the H/ emission (Section 3.1) and great progress has 
been made in observing and modeling the ionized gas 
(Section 3.3). However, the most significant step forward 
has probably been in our understanding of the physics 
of dust grains. With respect to the color variations 
mentioned above, PRONAOS has also shed light on the 
dust emissivity index /3: a significant inverse correlation 
between /? and dust temperature has been observed in 
the direction of star forming regions (Dupac et al. 2001). 
At the same time, analysis of the FIRAS data shows 
that, in the Galactic plane region, (3 is best-fitted by a 
value of 1 (Reach et al. 1995) rather than 2 as usually 
assumed. These are very important findings that alone 
provide suflacient motivation for revisiting the problem 
of decomposing the Galactic infrared emission. Further 
motivation is provided by the impending Planck mission 
(http//www.rssd.esa.int/ planck) which will obtain high- 
sensitivity, high-resolution (~ 5 to 30 arcmin) all-sky 
maps in the wavelength range ~ 400 /xm - 10 mm. Due to 
its unique instrumental performance, Planck will provide 
unprecedented insights on the radio/infrared emission 
of the diffuse ISM of our Galaxy. The present work was 
carried out as part of a larger effort to develop tools for 
the full exploitation of the Planck data. 

The paper is organized as follows. Section 2 describes 
the inversion technique, including a brief review of 
previous works. In Section 3 we provide details regarding 
the datasets used, with particular attention given to the 
ionized gas. Section 4 presents our main results while 
Section 5 provides conclusions and perspectives. 



2. The inversion technique 

Along the Galactic plane, the detected infrared emission is 
a blend of radiation arising from dust that is spread over 
a wide range of distances and Galactic radii, and subject 
to a wide range of physical conditions. The purpose of 
the present study is to decompose this integrated infrared 
emission into radial bins associated with each phase of the 
interstellar gas and determine the physical properties of 
each phase in each bin. This goal can be achieved by means 
of a so-called inversion method that employs kinematic 
distances to assign gas to radial bins and require that the 
gas distributions in the adopted bins have distinctly differ- 
ent spatial distributions. The inversion technique has been 
first applied to astrophysical data by Bloemen et al. (1986) 
to determine the Galactocentric distribution of gamma- 
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ray emissivity. Subsequently, Bloemen et al. (1990), Giard 
et al. (1994) and Sodroski et al. (1997) have used the same 
method to reconstruct the radial distribution of infrared 
emission. In the following, we will adopt the mathemati- 
cal formalism introduced by these authors to describe in 
detail this technique. 

Let Ix be the observed emission at wavelength A for a 
pixel j of a map. Then the modeled emission Im\ for this 
pixel can be written as: 



(1) 



i=l 



In the expression above: i denotes the intervals (or rings) 
of galactocentric radii Ri over which the decomposition is 
performed; chi {Ri, A), {Ri, A), {Ri, A) are the dust 
cmissivitics associated with the different phases of the gas 
in each ring; represents the H atom column density 
for neutral atomic hydrogen (H/) in the interval i?,; sim- 
ilarly, N^j^ is the H atom column density for molecular 
hydrogen (H2) and Nj^^^ is the H atom column density 
for ionized hydrogen (H//) in the considered ring. 

Note that eq. (1) does not require resolution of the 
kinematic distance ambiguity for material within the solar 
circle, since we assume that the dust properties vary only 
with Galactic radius. In addition, following Giard et al. 
(1994) and Sodroski et al. (1997) we do not assume a 
constant ratio between €Hi{Ri,^) and eH^{Ri,X), rather 
these parameters arc; allowed to vary in c;ac:li galactocentric 
ring. The dust emissivities can be determined by means 
of a least-square fit analysis, i.e., by minimization of the 
quantity: 



j=pixels 



{Im{j)^ 



I{j)x? 



(2) 



over the map. In the expression above, Uj {uj = Uj{X)) 
represents the noise per pixel for the I\ map (see Section 
4.1 for details). 

Column densities are computed, for every ring and for each 
phase of the gas, by applying a circular rotation model for 
the Galaxy. The rotational velocity 0{R) of a gas clement 
in differential motion with respect to the Galactic center 
is given by: 



e{R) 



R 
Ro 



^0 + 



Vlsr 
sin I cos 6 



[km/s] (3) 



with R the galactocentric distance of the considered gas 
element, Rq the galactocentric distance of the Sun (we 
take Rq = 8.5 kpc), 9o the rotational velocity of the Sun 
(we assume — 220 km/s), Vlsr the linear velocity along 
the line of sight in the Local Standard of Rest and I and 
b the Galactic coordinates for a given direction. We adopt 



the rotation curve e{R) by Fich, Blitz & Stark (1989) for 
which the galactocentric distance can be expressed as: 



R = 



221.641 X Ro 



60 -h 0.44286 X Rq 



Vls 



[kpc] (4) 



sin / cos h , 



Therefore, given a specified ring Ri, the corresponding 
range in linear velocities is: 



LSR 



221.641 X 



Ro 

Ri 



[km/s] 



(5) 



sinlcos;(6'o + 0.44286 x Ro) 



For each given ring, iV|^^, Njj-^ and Njj-^^ are then com- 



puted by integration over the velocity channels selected 
as above. Further details, specific to each gas phase, will 
be given in the following section. We note that the use 
of a different rotation model would not significantly affect 
the results of the analysis. For the ring decomposition, we 
have chosen the following intervals: 



Ri = [[0.1,4], [4,5.6], [5.6,7.2], [7.2,8.9], 
[8.9, 14], [14, 17]] 



[kpc] (6) 



The above decomposition has been obtained through 
a process of optimization of the number and type of 
Galactocentric intervals aimed at reducing the cross- 
correlation between neighboring rings and therefore lead- 
ing to a well-defined linear inversion problem. This goal 
has been achieved by starting with a fine decomposition 
(roughly 200 bins in Galactic radius) and then gradiially 
combining correlated rings, guided by the eigenvalues of 
the cross-correlation matrix for different components. 
All the maps used for the inversion have been projected 
into HEALPix format (Gorski et al. 1999) and degraded 
to 1° resolution. The final maps (input IR maps as well 
as H/, II2 and H// galactocentric-annuli maps) arc at the 
HEALPix resolution (i.e. nside) 256 which roughly corre- 
sponds to a pixel linear size of 13.7'. 

3. The data base 

3.1. Dust 

To trace dust emission we use a large data base ranging 
from sub-millimeter to centimeter wavelengths. At short 
wavelengths, i.e. 60 and 100 /xm, we use the last generation 
of IRAS plates, the so-called IRIS {Improved Reprocessing 
of the IRAS Survey) maps (Miville-Deschenes & Lagache, 
2005). These maps have been obtained by reprocessing 
the original ISSA plates in a manner that allowed the 
correction for calibration, zero level and striping prob- 
lems. The IRIS maps (which are at an angular resolution 
of 4-arcmin at 60 /im and 4.3-arcmin at 100 /im) also 
have a better zodiacal light subtraction. At 140 and 240 
/xm, we use the DIRBE Zodi- Subtracted Mission Average 
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Maps. DIRBE is a photometer on-board the COBE satel- 
lite with a resolution of 40-arcmin. To probe the sub- 
millimetcr wavelengths, we make use of Archcops data. 
Archeops is a ballon-borne experiment (Benoit et al. 
2002) which has detected Cosmic Microwave Background 
(CMB) anisotropies at the average angular resolution of 
8-arcmin in 4 frequency bands, 143, 217, 353 and 545 GHz. 
It is the only experiment of this kind which has provided 
information on the Galactic plane emission at millimeter 
wavelengths over a large fraction of the sky. Three flights 
have been performed: a first one (T99) in July 1999 from 
Trapani (Italy) followed by a second one (KSl) in January 
2001 from Kiruna (Sweden) and by a third one (KS3) in 
February 2002 again from Kiruna. For the present anal- 
ysis, we use the T99-KS3 combined data at 143 and 353 
GHz. These combined maps provide better sky coverage 
than individual maps. In particular, the T99-KS3 data set 
covers the longitude interval 30 < 1 < 180 degrees. At 
centimeter wavelengths, we use the WMAP all-sky maps. 
WMAP is a satellite dedicated to the observations of CMB 
anisotropies with a frequency coverage in the range 23 
to 94 GHz and an angular resolution increasing with fre- 
quency from 52.8 to 13.2 arcmin. 

3.2. Atomic hydrogen 

For the atomic hydrogen, we use the newly-released 
LAB (Leiden/ Argentine/Bonn) survey (Kalberla et al. 
2005). This data set has been obtained by merging the 
Leidcn/Dwingeloo Survey (hereafter LDS, Hartmann 
& Burton, 1997) of the sky north oi 5 = -30° with 
the Institute Argcntino dc Radioastronomia Survey 
(hereafter lAR, Arnal ct al. 2000; Bakaka et al. 2005) 
which covers declinations south of (5 = -25°. The angular 
resolution of the final data cube is 0.6° while the velocity 
resolution is 1.3 km s~^ with a LSR velocity range from 
-450 km s~^ to -1-400 km The lAR survey was 

specifically designed to complement the LDS survey of 
the northern sky. Although the observational parameters 
of the two surveys are quite similar, some reprocessing 
has been necessary to combine them. Both data sets 
have been corrected for stray light radiation. In addition, 
the existing LDS cube has been reprocessed to correct 
for reflected ground radiation. Owing mainly to limited 
computing power, the original processing was not able 
to apply this correction accurately. Consistency checks 
between the two surveys have also been made in the 
region of overlay (-30° < 5 < -25°), and they appear 
to agree quite well (Bajaja et al. 2005). To merge the 
two data sets into a single cube a spline interpolation 
followed by a spatial regridding has been performed. Such 
regridding has slightly degraded the original resolution of 
the two surveys to 37.5'. The rms brightness-temperature 
noise of the merged database is 0.07-0.09 K. Residual 
errors due to defects in the correction for stray radiation 
are mainly below 20-40 mK. 



This data base is used to compute, for every ring Ri 
and every given line of sight k, the quantity {Nhi)\ ac- 
cording to (Kerr 1968): 

{N\j^)k = 1.82 X 10-'T« X (7) 

In the expression above, is the spin temperature 
(Ts=145 K) while T^, in K, is the H/ observed bright- 
ness temperature. The integral is taken over the velocity 
range defined in eq. (5). Units are 10^° atoms cm~^. It is 
important to emphasize that different assumed values of 
Tg in the range 120-160 K would change computed column 
densities by less than a few percent (Sodroski et al. 1997). 

Fig. 1 shows spatial maps of HI integrated over each 
of the radial bins. The emission has the widest latitude 
extent in the 7.2 < R < 8.9 kpc bin - which includes the 
Sun - and becomes systematically narrower in the smaller 
and larger radial bins. A feature of note in the 7.2 < R 
< 8.9 kpc bin, also seen in the corresponding H2 map in 
Fig. 2, is the lack of emission in the range -140 < 1 < -100 
deg. This gap may be a consequence of the structure of 
the Gould Belt and the corresponding Lindblad expanding 
ring. In the general direction of the gap, the Gould Belt is 
quite far from us, far enough that its gas would be shifted 
into the R=8.9-14 kpc bin (see, e.g., Fig. 5 in Perrot & 
Grenier 2003). In the highest radial bin (14 < R < 17 
kpc) one sees clearly the well-known warping of the outer 
Galactic disk. 

3.3. Molecular hydrogen 

Templates of the radial distribution of H2 column densities 
have been constructed from CO data. In particular, we use 
the Galactic CO survey of Dame et al. (2001) except in 
a region covering parts of the first and fourth quadrants 
( 0° < 1 < 60° and 273° < 1 < 360°, -5° < b < +4°) 
where the NANTEN survey (Fukui et al. 1999) is used. 
The NANTEN data set was obtained with the Nagoya 
University 4-m radio telescope located at Las Campanas 
in Chile. The final survey is a composite of 1,100,000 spec- 
tra at 4 spacing. The velocity resolution is 0.65 km/s al- 
though the version of the data we worked with has been 
slightly degraded to a resolution of 1 km/s. The veloc- 
ity range is -250 to -1-250 km/s. As for the Dame et al. 
(2001) data set, this has been obtained, although at dif- 
ferent times, with two very similar 1.2-m telescopes, one 
in Cambridge, Massachusetts and the other at the Cerro 
Tololo Inter- American Observatory, Chile. The final com- 
posite CO data set has been constructed by interpolating 
all component surveys onto a uniform 1-b-v cube with an- 
gular spacing of 0.125° and a velocity spacing of 1.3 km/s. 
In addition, bad channels and single missing spectra were 
filled in by interpolation. The new CO survey has 16 times 
more spectra than that of Dame et al. (1987), 3.4 times 
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Longitude [deg] 

Fig. 1. Ring 1 to 6 (from top to bottom) for H/ column densities. Levels of the contours are: 3, 30, 70, 110 and 150 
10^° atoms cm~^. 



higher angular resolution and up to 10 times higher sensi- 
tivity per unit solid angle. The proton column density in 
molecular hydrogen for a given ring Ri and line of sight k 
is then obtained (Lebrun et al. 1983) as: 

(A^l^Jfc = 2(2.8 X 10^0) /" Tcoik,v)dv (8) 

Units are in 10^" atoms cm~^. In the expression above, 
following Giard et al. (1994), we assume a Galaxy- wide 
constant value of 2.8 x 10^° mol cm~^ km~^ s for the 



ratio of H2 column for the ratio of H2 column density to 
velocity-integrated CO intensity (Bloemen et al. 1986). In 
addition, the integral is computed over the velocity range 
defined in eq. (5). 

Fig. 2 shows the spatial maps of CO integrated over 
the same radial bins as the H/ in Fig. 1. The most striking 
difference between the two figures is the lack of CO emis- 
sion in the highest radial bin, demonstrating the much 
larger radial extent of the Galactic H/. Evident in the CO 
local bin (7.2 < R < 8.9 kpc) are well known molecular 
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Fig. 2. Ring 1 to 6 (from top to bottom) for CO column densities. Levels of the contours are: 20, 40, 50, 90, 130 and 
170 10^0 atoms cm'^. 



clouds such as the Cepheus flare (/ ~ 100 deg), the Aquila 
Rift (Z ~ 20 deg) and the Chamaeleon (/ ~ -60 deg). The 
Orion complex appears mainly in the 8.9 < R < 14 kpc 
bin. 

3 A. Ionized hydrogen 
3.4.1. Previous works 

Before illustrating how the ionized gas has been handled 
in our analysis, we briefly review the approaches adopted 



by previous authors. Over the years, the contribution of 
ionized gas to the observed IR dust emission has been the 
subject of considerable debate. According to Mezger et 
al. (1982) and Cox ct al. (1986) 60% to 80% of the total 
infrared luminosity of the Galactic disk is contributed by 
dust associated with ionized gas. Sodroski et al. (1989) 
in their early analysis attempted to model this contri- 
bution by means of the 5 GHz radio continuum data of 
Haynes et al. (1978). To account for the non-thermal syn- 
chrotron component at this frequency, they estimated the 
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Fig. 3. Ring 1 to 6 (from top to bottom) for Hjj'^^ff'^^'^ column density. Levels of the contours are: 
5,20,40,80,120,160,200,300 10^° atoms cm'^. 



thermal/non-thermal contributions at 1.4 GHz and then 
extrapolate that contribution to 5 GHz by adopted con- 
stant spectral indices, i.e. -0.1 for the thermal component 
and -0.7 for the non-thermal on43- They concluded that 
less than 20% of the 5 GHz emission is of non-thermal 
origin. However, this method has a major problem. As 
shown by Paladini et al. (2005), a careful separation of the 
thermal and non-thermal components has to take into ac- 
count the spatial variation of the synchrotron spectral in- 

^ Quoted spectral indices assume S„ oc i/". 



dex. When such variation is considered, one finds that the 
non-thermal emission contributes only ~ 20% at |6| 1°. 
At higher latitudes, the synchrotron component increases 
significantly. In addition, it is shown that discrete H// re- 
gions contribute about 10% of the total thermal emission. 
A slightly different approach is adopted by Sodroski et al. 
(1997). In their paper, the 5 GHz data are complemented 
with 19 GHz low-resolution (180 ) data from Boughn et al. 
(1991). As shown by Finkbeiner et al. (2004), the emission 
outside of H// regions for v >\A GHz is likely dominated 
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by spinning dust. Therefore the use of data at 19 GHz as 
a tracer of the ionized gas turns out to be rather problem- 
atic. On the contrary, Blocnicn ct al. (1990) cxchidc the 
ionized gas phase from their model based on the work by 
Cox and Mezger (1988), who concluded that dust associ- 
ated with ionized gas contributes only 20% to the overall 
IR luminosity. As a consequence, bright H// regions are 
masked on the basis of their 60/im/lOO/im colors and ex- 
cluded from the fit. However, as noted by Giard ct al. 
(1994), the large discrepancy between the IR observations 
and their best-fit model intensities is likely due to the ex- 
clusion of the ionized gas component from their analysis. 

3.4.2. This work 

The best way to trace the ionized gas and its spatial 
distribution in the Milky Way would be a Galactic plane 
survey of radio recombination lines (RRLs). Since such a 
survey is not yet available, we have to resort to alternative 
methods. Taylor & Cordes (1993) have proposed a model 
which is based on observations of pulsar dispersion and 
scattering measures (respectively, DM and SM). The 
model has been incorporated and improved by Cordes & 
Lazio (2002, 2003, hereafter NE2001). The NE2001 model 
consists of five components: a large scale distribution, the 
Galactic Center, the local ISM, individual clumps and 
individual voids. The large scale distribution is, in turn, 
contributed by two axisymmetric components, i.e. a thin 
and a thick disk, and by the spiral arms. For details about 
the functional form of these components, we refer the 
reader to Cordes & Lazio (2002). By using the NE2001 
model, we can construct rings (see Fig. 3) of electron 
column densities for the diffuse ionized gas following 
the same guidelines as for the H/ and H2. In order to 
check the reliability of this model, it would be highly 
desirable to tie up the computed column densities with 
observations. This corresponds in practice to converting 
column densities (DM) into emission measures (EM). 
Although some information exists on the relation between 
these two quantities (sec, for example, Mitra, Berkhuijsen, 
MuUer (2003)), it is still too uncertain to allow a realistic 
conversion and therefore to validate and confidently use 
the model. 

Another way to trace the ionized gas is through its 
free- free emission. With the release of WMAP data, all-sky 
templates of the Galactic foregrounds have become avail- 
able (Bennett et al., 2003b). These include synchrotron, 
free-free and dust emission maps. Such maps have been 
obtained with the Maximum Entropy Method (MEM). 
With this technique, prior templates for each component 
of Galactic emission are considered and it is assumed that 
the observed signal at each frequency is contributed by 
four terms, namely, CMB, synchrotron, free-free and dust. 
For free-free emission (i.e. thermal brcmstrahlung emis- 
sion produced by ionized gas), the Hcc map of Finkbeiner 
(2003) is used as a template. We note that Ha maps can- 



not be used directly to trace the ionized gas given that 
they suffer heavily from obscuration due to dust absorp- 
tion at low latitudes (|fe| < 5 deg). As reported in Bennett 
et al. 2003b, all the uncertanties associated with the Ha 
map (i.e., errors due to correction for dust as well as to 
conversion from Ha to free- free, etc.,) combine together 
and do not allow one to put tight constraints on the recov- 
ered free-free map. Despite these drawbacks and the fact 
that the WMAP free-free map traces continuum emission 
and, as such, does not allow the recovery of the spatial 
location of the emitting gas, we have decided to make use 
of it for our analysis, convinced that it represents the best 
available option. The WMAP team has released free-free 
maps at each WMAP band. We have used the one in the 
Ka band (i.e. 33 GHz) as a compromise between spec- 
tral dependence and angular resolution. Such a map, pro- 
vided in units of thermodynamic mK, has been converted 
into imits of 10^" electrons cm~^ by applying eq. (8) in 
Sodroski et al. (1989) with a value of 10 cm~^ assumed for 
Ueff- We have compared the performance of the WMAP 
free-free map with respect to the NE2001 model. The re- 
sult, in terms of reduced x^, i.e. x^, for a sample of wave- 
lengths is shown in Table 1. 





60 /um 


100 urn 


140 urn 


240 urn 


WMAP free-free 


1.087 


1.046 


1.037 


1.033 


NE2001 


1.141 


1.078 


1.056 


1.043 



Table 1. comparison for WMAP free- free map and 
NE2001 electron density model. 



Clearly, the fit to the observed omission gets worse when 
the NE2001 model is applied. This is expected since the 
NE2001 model mostly accounts for the very large-scale 
structure while, as noted by Sodroski et al. (1989), the 
bulk of the IR emission consists of a distribution of 
prominent, discrete sources superposed on a continuous 
background. 

The WMAP free-free map is characterized by a rather 
low angular resolution (1°). This implies that most of the 
sources (H// regions) which contribute to the Galactic 
thermal emission are highly diluted since their average 
angular size is ~ 6 . Although the diffuse component ac- 
counts for nearly all of the ionized gas in the Galactic disk 
(its mass is a factor of 20 larger than that of classical H/j 
regions (Reynolds 1991)), we have included the contribu- 
tion of discrete sources for completeness. For this purpose, 
we have made use of the Paladini et al. (2003) catalog 
(hereafter Paper I). This catalog contains 1442 Galactic 
H// regions for which angular diameter and flux densities 
at a reference frequency (2.7 GHz) are given. Radio re- 
combination lines have been observed for ~ 800 of these. 
This kinematic information, combined with the Fich, Blitz 
& Stark (1989) rotation curve, has been used in Paladini, 
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Davies & DeZotti (2004) (hereafter Paper II) to compute 
galactocentric and solar distances. Only sources (575) for 
which the observed velocity is |Kfcs| > 10 km/s have been 
considered, to minimize uncertainties in the derived dis- 
tances due to peculiar motions. For 281 H// regions lo- 
cated inside the solar circle, the distance ambiguity has 
been resolved by applying one of these methods: 

1. comparison with auxiliary (optical or absortion) data 
when these are available 

2. use of a statistically significant luminosity-physical di- 
ameter correlatioro 

Following these guidelines, we obtain 550 H// regions for 
which the solar distance is uniquely determined. In Paper I 
we have also investigated the completeness of the data set. 
Such a study is complicated by the fact that the catalog 
is constructed by combining data obtained with different 
instruments, observing at different frequencies and angu- 
lar resolution. Based on simple geometric arguments, we 
estimate that the catalog is complete, at 2.7 GHz, for flux 
densities > 7 Jy (see Fig. 1 of Paper I). For fainter sources, 
the catalog completeness is mostly limited by confusion. 
For each source, we derive the electron density according 
to (Schraml & Mezger, 1969): 

ne = 98.152 X T^-^^^ S"-' D""'^ 0-^-^ [cm'^] (9) 

being v the frequency in GHz, the electron tempera- 
ture in K, S the flux density in Jy, D the solar distance 
in kpc and 9 the angular size in arcmin. From we then 
compute the electron column density simply by multiply- 
ing for the linear diameter expressed in pc. For most of 
the sources, an electron temperature has been computed 
from the observed recombination line. However, when Tg 
is not available, we estimate it by applying the empirical 
relation: 

re(i?) = (4166 ± 124) + (314 ±20)i? [K] (10) 

as derived in Paper II, with R, the galactocentric radius, 
expressed in kpc. As can be seen in Fig. 4 of Paper II, the 
Te — R correlation is affected by a large scatter. However, 
eq. (2) has been used only for 20% of the sources. 

4. Derived physical parameters 

4.1. The iterative inversion method 

In matrix form, the minimization of eq. (2) is equivalent 
to solving the linear system: 

{A'nA)ex = A'nh (11) 
where A is a (3n 4- 1) x matrix of the type: 

rli rli n2 rl2 n.ii J^ii J^ii 

Ul Hi U2 rl2 rill ri i j Uu 

^ See eq. (6) in Paper II 



n being the number of rings in the decomposition and 
A^ the number of pixels. In the expression above, N^^''^ 
denotes the column density in pixel j and ring i for the 
source (or compact) component of the ionized gas, while 
Nj^^^ is the corresponding column density for the diffuse 
component. H has dimension N x N and is given by: 

° 
V ... -i^ / 

with aj the pixel noise of the input map. In eq. (11), 
the correlation matrix A*nA is symmetric and positive- 
defined which allows the linear system to be solved by 
means of a Cholesky inversion method (see, for instance, 
Golub & van Loan, 1989). The major difficulty is repre- 
sented by a correct evaluation of the noise which is largely 
dominated by the uncertainty associated with the model. 
To circumvent such a problem, we have applied an iter- 
ative method. This consists of making an initial assump- 
tion about the (instrumental) noise, computing the corre- 
sponding solution and then using the noise derived from 
this partial solution to reiterate the process. A true esti- 
mate of the real noise is obtained when, by repeating this 
process several times, the system converges. Convergence 
is, on average, reached within 10 iterations. In the first 
iteration the (intrumental) noise per pixel is assumed to 
be given by: 

+ 5%/a(j) (12) 

where ai„i represents the average noise of the Ia map at 
high latitudes and lA(j) is the value of the map in pixel j. 
In this expression, the term proportional to the observed 
signal is very important. The sky is in fact highly inhomo- 
geneous and variations in signal in the range 1 to 10,000 
(in brightness units) are observed when high latitutes are 
compared to low ones. It is then very unlikely that the 
difference between the true and the measured signal is 
only the same additive noise as it is at high Galactic lat- 
itudes. A proportional term taking into account the long 
term drift of the instrument gain needs to be considered. 
Nonetheless, it is important to emphasize that such defini- 
tion of the noise is still rather crude and purely operational 
given the iterative method which has been adopted. 
Emissivities coefficients obtained by following the proce- 
dure now described are shown in Table 2. 
Uncertainties quoted in the table have been computed by 
applying the bootstrap method described in Giard et al. 
(1994). 

4.2. Results and discussion 

Fig. 4 and Fig. 5 show longitude profiles of the ob- 
served and modeled intensities as well as the respective 
contributions to the calculated profiles from H/, H2, 
diffuse Rj/'ff^'^^ and H//'^"'"'"''^*. Fig. 6 illustrates the 
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Fig. 4. Longitude profiles for IRAS 60 and 100 /im (top panels) and DIRBE 140 and 240 fim (bottom panels). Black 
line denotes observed emission while light green line corresponds to the fitted model. In addition, blue line is for H/, 
dark green line for H2, red line for H//*-''^"*'^ and magenta line for HjjCompact ^ Profiles are obtained by averaging over 
the shown latitude intervals. The apparent noise is not instrumental: along the Galactic plane, the S/N is very high 
and error bars representing the instrumental noise would be of the same size of the plotting line. 
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Archeops 353 GHz : 



180 160140120100 80 60 40 20 -20-40-60-80-10(^12O14^16i}180 
Goladic longitude (degree) 




Archeops 143 GHz 



180 160140120100 80 60 40 20 -2O-4O-6O-8O-lOO120-14^16i}18O 
Goloctic longitude (degree) 




Fig. 5. Longitude profiles for Archeops 850 and 2096 (top panels) and WMAP W and Q bands (bottom panels). 
Archeops data base covers only the longitude interval 30° < / < 180°. Black line denotes observed emission while 
light green line corresponds to the fitted model. In addition, blue line is for H/, dark green line for H2, red line for 
jljjdiffuse g^^^ magenta hue for JijjCompact^ Profiles are obtained by averaging over the shown latitude intervals. The 
same remark as in Fig. 4 on instrumental noise applies here. 
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Fig. 7. Galactoccntric distributions of the infrared emis- 
sivities for dust associated with H/ (top left), H2 (top 
right) and H//'^°'"*^'^* (bottom center). Archeops results 
at 850 /im have been obtained for the longitude interval 



corresponding latitude profiles. The fitted region in the 
longitude profiles, at each frequency, covers the latitude 
range \b\ < 5°, due to the limited sky coverage of the CO 
observations. We also excluded the longitude range \l\ < 
5° and the corresponding region in the anticenter due 
to the lack of kinematic resolution in these directions. 
Although we could have corrected for this effect by, for 
instance, constraining the radial velocity interval, follow- 
ing the method implemented by Bcrtsch ct al. (1993), we 
decided to ignore such a problem in the current analysis 
and to investigate the Galactic center region in detail in 
a forthcoming paper. 

The longitude profiles show that the model fits are 

satisfactorily good at all wavelengths with the exception 
of a few peaks at 60 /um (which are not correctly repro- 
duced) and some overshooting in the WMAP bands. The 
peaks are likely H/j regions and the model fails in this 
case because the H// region catalog of Paladini et al. is 
complete only for fairly bright sources (see Section 3.4.2). 
In addition, as discussed in Section 3.4.2, the angular res- 
olution of the WMAP free-free map does not allow one to 
resolve discrete sources. The overshooting testifies to an 
excess of emission and this can be explained by the fact 
that WMAP channels are mostly dominated by thermal 
bremstrahlung instead of dust emission (see also Fig. 8 
and related discussion). 

Similarly, the latitude profiles show that the model re- 
produces quite well the observed emission. Moreover, one 
can see that the emission between 60 and 2096 /im is dom- 
inated by dust associated with H/ although at 60 /im this 
is largely contributed by dust associated with the ionized 
gas. Remarkably, the H// region contribution is clearly 
visible in the IRAS 60- /im profile but almost disappears 
in the corresponding plot for WMAP at 41 GHz despite 
the fact that tlic^ emission, at this wavelength, is mostly 
thermal radiation from the ionized gas. As discussed pre- 
viously, this is an effect due to the angular resolution: 
IRAS resolving power matches the average angular size of 
cataloged H/j regions while WMAP resolution is 10 times 
lower. 

Fig. 7 shows the Galactoccntric distribution of dust 
emissivities associated with H/, H2 and H// regions. 
The case of diffuse H// is not plotted since for that 
component we do not have any 3d-spatial information. 
The emissivities for dust associated with H/ decreases 
strongly with increasing Galactoccntric distance at all 
wavelengths with the exception of 850 /xm which shows 
a less pronounced trend. Such behaviour appears both 
inside (R < Rq) and outside (R > Rq) the solar circle. On 
the contrary, the radial distribution for the emissivities 
of dust associated with H2 shows a bump, at all A, at R 
~ 5-6 kpc. Similarly, the distribution for dust associated 
with H// regions presents a first pronounced bump at R ~ 
5 kpc and a second at R ~ 8 kpc. In the case of H/, if we 
combine this information with the fact that the intensity 
of the ISRF scales with Galactoccntric distance in a 
comparable way (Mathis et al. 1983), we can conclude 
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that dust mixed with atomic hydrogen is mostly heated 
by the general radiation field, in agreement with previous 
findings by Bloemen et al. (1990) and Sodroski et al. 
(1997). The behaviour of dust associated with molecular 
hydrogen suggests instead a correlation with the intense 
star formation activity taking place in the molecular ring, 
i.e. dust associated with H2 appears to be heated also 
(and perhaps by a large fraction) by massive O and B 
stars still embedded in the parent molecular clouds rather 
by only the general radiation field . Sodroski et al. (1997) 
reached a similar conclusion while Blomen et al. (1990) 
did not find any significant correlation between the radial 
distribution of dust emissivity associated with H2 and 
the molecular ring. As for H// regions, the derived radial 
distribution of emissivities is fully consistent with the 
radial distribution of sources found in Paladini et al. 
(2004; see their Fig. 3). 

The recovered emissivity coefficients allow one to com- 
pute the 60 /im/100 /im ratio in each Galactocentric bin, 
as shown in Table 3. The ratio is fairly constant for dust 
mixed with molecular gas with the exception of the first 
ring which is however characterized by a significant er- 
ror. In the case of dust mixed with atomic gas though, 
we find a significant decreasing trend with Galactocentric 
radius although, within the solar circle (R < 8.5 kpc), 
the gradient appears to be quite small. This second re- 
sult partly contradicts Bloemen et al.'s finding that the 
60 /im/100 /im ratio is nearly constant up to R ^ 17 kpc. 
Nevertheless, when we take the mean values, we obtain 
60 /xm/lOO /J,m ratios of ^ 0.22 and 0.23 for, respectively, 
dust associated with atomic and molecular ga^ in agree- 
ment with Bloemen et al. (who report average values of 
~ 0.27 and 0.20) and also with typical values at high lati- 
tudes (for instance, Boulanger & Perault (1988) derive an 
average value of 0.21±0.02 at |6| < 10 deg). 
The derived IR emission spectra for dust associated with 
each phase of the gas are shown in Fig. 8. We have fitted 
the FIR luminosity per H mass, Airvei, (in units of 10"'^^ 
W (H ATOM)-i), with a modified-blackbody, i.e.: 

47rzye^ = 47rAeA ^ AirXaBTiX) f ^ j (13) 

In this expression, a represents the dust cross-section per 
H atom (in units of / H atom). 

We have assumed a spectral emissivity index fi of 1.5. 
Such a value represents a compromise between the results 
obtained by Dupac et al. (2003) for M17, i.e. 1.3 < /3 < 1.8 
in the temperature range 35-15 K, and the fact that /3=1.7 
is usually adopted in the Solar neighborhood (Lagache et 
al. 2003). This flat spectral index is indicative of the exis- 
tence of a millimeter excess with respect to the (3 = 2 be- 
havior expected from several grain models. Such an excess 
was observed for the first time in analyzing FIRAS data at 
A > 500 nm (Wright et al. 1991). A similar finding is also 

^ The latter has been computed without the value for the 
first ring. 



supported by observation of external galaxies as shown by 
GaUiano et al. (2005). Finkbeiner et al. (1999) attribute 
this excess to a very low temperature (9 K) emission com- 
ponent in the diffuse ISM. Alternatively, the excess could 
be attributed to a temperature dependence of the dust 
submm emissivity spectral index (Meny et al. 2006) 

By fitting the spectra with eq. (13), we derive esti- 
mates of the dust temperatures for every gas phase and for 
each radial interval (see Table 4). The following average 
temperatures have been obtained: T(;(H7) = 19.86±2.84 
K, Td(H2)=19.23±2.ldl K, Td(H//''*-''/"''^)=26.7±0.1 K 
and Td(H/7™™P"'=*)=27.87±2.89 K. 

Clearly, the WMAP points (with the exception of W 
band) cannot be fitted by means of a modified black-body. 
Likely, this wavelength range is dominated by a mixture of 
thermal bremstrahlung (free- free), synchrotron emission 
and possibly (see Section 3.4) spinning dust rather than 
thermal radiation from stable dust in the diffuse ISM. This 
combination of radiative processes is not easy to model 
since the 3d-templates are not available and the synchro- 
ton spectral index varies spatially. Therefore, we fit the 
WMAP V to K bands with a power-law of the form a 

where a is a constant of normalization and /3 the 
spectral index for such a mixture of emission mechanisms. 
The power-law fit appears to represent satisfactorily the 
emissivity coefficients derived through the inversion tech- 
nique. The best-fit a and P values are shown in Table 4. 

An additional interesting result of the fitting pro- 
cedure is represented by the analysis of the Archeops 
point at 2096 /im as well as by the WMAP W-band 
value. In the outer Galaxy (Fig. 8, Ring 5 and 6) these 
points show a tendancy to lie above the fitted modified 
black body. This indicates that the galactic millimeter 
excess emission might be comparatively stronger towards 
the outer regions of the Galaxy. A similar finding was 
reported by Bourdin et al. (2002) also based on the use 
of the Archeops data. It is important to emphasize here 
that this result cannot be an artifact of our data analysis 
techniques. Such a spurious effect might be a matter of 
concern for the Archeops data set: the maps produced 
by the Trapani and Kiruna flights (Delabrouille, J. & 
Filliatre, P., 2004) have in fact been processed in a slightly 
different way, the Trapani map being signiflcantly more 
filtered than the Kiruna map to correct for systematics. 
However, since such filtering affects the larger rather 
than the smaller scales, it would tend to decrease rather 
than increase the observed excess. At the same time, 
we point out that this result has been obtained by fix- 
ing /3 to 1.5. A /3 of 2 would enhance the excess even more. 

The emissivity coefficients can also be combined with 
information on the mass distribution of the gas to compute 
FIR luminosities. In particular, if we denote with i a given 



^ This value is computed by averaging ring 2 to 5. Ring 1 
and 6 have been excluded given the significant uncertainty. 
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phase of the gas, the total FIR luminosity for this phase 
is obtained according to: 

LFiR,i = Mi / 47rAei,AdA [IO^Lq] (14) 

where is the mass of the gas (expressed in 10^ Mq), 
Attvcx is the FIR luminosity per H mass (in L0/M0 units) 
and the integral is taken over the interval [Xmin,^max] 
= [60 /im, 2096 fim\. In writing cq. (14) we use again 
the fact that Ani'ei^i, = 47rAei,A- As well as for the total 
Galaxy, FIR luminosities have been computed, for each 
phase of the gas, in each ring. In this case, the quantity 
Mi represents the mass of the gas in the considered 
Galactocentric interval. Masses for the atomic and 
molecular hydrogen have been derived by integrating, 
for each ring, the surface densities quoted by Binney & 
Merrifield (1998). As far as the ionized gas is concerned, 
we have made direct use of the estimate of the total mass 
reported by Westerhout (1958) for an effective electron 
density ne=10 cm~^. Remarkably, the value provided by 
Westerhout is a factor ^ 40 lower than the one given in 
Ferriere (2001). The two values have been obtained in a 
different way: the Westerhout's estimate is derived from 
free-free emission observations which arc proportional to 
the emission measure (EM) while the Ferriere's value is 
based on the Cordes et al. (1991) model which provides 
a direct measure of the electron column density (DM) 
(Katia Ferriere, private communication). The fact that 
our Galaxy contains several thermal components which 
are characterized by different clumping factors, i.e. 
different relations between DM and EM, explains the 
discrepancy in the two estimates. 
The resulting FIR luminosities are given in Table 5. 

Finally, we have computed the contribution of each 
component to the global SED. Fig. 9 shows that the dom- 
inant contribution at all wavelengths > 60 /um is due to 
dust associated with atomic gas. However, at A ~ 60 fim 
the emission is largely contributed by dust in H77 regions. 



5. Conclusions 

By using an inversion technique, we have derived the ra- 
dial distribution of dust properties in our Galaxy in the 
wavelength range 60 to 2096 /im. In particular, we have 
obtained emissivity coefficients for dust associated with 
atomic, molecular and ionized gas. By assuming a modi- 
fied black-body with /3=1.5, we estimate that average tem- 
peratures for dust uniformly mixed with H/, H2 and H// 
are of order of 19.8, 19.2 and 26.7 K respectively. It is im- 
portant to emphasize that the derivation, through inver- 
sion techniques, of dust emissivity coefficients associated 
with the molecular and ionized phases may be affected by 
a mutual contamination due to the spatial correlation that 
exists between these phases. In turn, such contamination 
would affect the derived physical quantities such as tem- 
perature. As noted by Sodroski et al. (1997), only the use 



of RRLs to trace the ionized gas phase would allow one to 
circumvent this problem. From the radial distribution of 
the emissivity coefficients associated with each gas phase, 
we find evidence that dust associated with the H/ is heated 
by the global ISRF. On the contrary, dust in molecular 
clouds appears to be heated in a significant way by young 
massive stars still embedded in their parent clouds. Our 
best-fitted spectra show an excess at Archeops 2096 /xm 
and WMAP 94 GHz. Such an excess, also found in FIRAS 
data, can be interpreted as due to a very cold dust compo- 
nent in the diffuse ISM or to a temperature dependence 
of the spectral emissivity index at submm wavelengths. 
In addition. The long-wavelength (A > 4900 fim) ) part 
of the spectra for the H/ and H2 associated dust compo- 
nents appears to be dominated by a mixture of thermal 
bremstrahlung, synchrotron emission and likely spinning 
dust. The possibility of a significant contribution due to 
spinning dust will be analyzed in a forthcoming paper by 
means of auxiliary data at even longer wavelengths (A > 
13000 /im). As shown by recent studies (i.e. Watson et al., 
2005), the peak of anomalous emission is expected to be 
at 15000 Atm. 

The derived emissivity coefficients also allow us to com- 
pute FIR luminosities. For these, the dominant contribu- 
tion (~ 80%) results to be provided by dust associated 
with atomic gas. 
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Table 2. Conversion factors to Aniye^ are, in units of 10"^^ W (H atom)"^: 6.28, 3.77, 2.68, 1.57, 0.44, 0.18 at 60, 
100, 140, 240, 850 and 2096 fim respectively. Quoted errors are derived as described in Section 4. 



Gas phase Ring eeonm/eiooum 
(kp£) 

H/ 0.1-4 0.28±0.02 

4-5.6 0.28+0.01 
5.6-7.2 0.25+0.01 

7.2-8.9 0.23±0.002 

8.9-14 0.16±0.003 
14-17 0.14+0.06 

Hz 0.1-4 0.7±0.51 
4-5.6 0.27±0.02 
5.6-7.2 0.22+0.02 
7.2-8.9 0.22+0.01 
8.9-14 0.2±0.01 
14-17 



Table 3. 60 over 100 /um ratio for dust associated with H/ and H2. 
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Fig. 6. Latitude profiles (from top left, clockwise) for: IRAS 60 and 100 ^m, DIRBE 140 and 240 /^m, Archeops 850 
and 2096 fim, WMAP W and Q bands. Black line denotes observed emission while light green line corresponds to 
the fitted model. In addition, blue line is for H/, dark green line for H2, red line for Hjj^-^ff'^^'' and magenta line for 

^compact 
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Fig. 8. Derived emission spectra for dust associated with H/ (top left), H2 (top right), H//^'""^"'^* (bottom left) and 
]^jjd.iffuse (ijottom right). Dashed lines illustrates best-fitted modified black-body and power-law functions. Symbol 
notation is as in Fig. 7: the star corresponds to IRAS 60 ^m, the square to IRAS 100 /im, the triangle to DIRBE 140 
/xm, the asterisk to DIRBE 240 fim, the filled circle to Archeops 850 /xm and the black diamond to WMAP W-band. 
In addition, red diamonds denote WMAP K, Ka, Q and V bands. 
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Table 4. Derived temperatures for dust associated with each gas phase at increasing radial distance from the Galactic 
center. 
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Table 5. Emissivities per H mass (L'FiR,i/^i) are computed from values of quoted in Table 2 by using the following 
conversion factors: 1.92 at 60 /xm, 1.15 at 100 /xm, 0.82 at 140 /xm, 0.48 at 240 /xm, 0.13 at 850 /xm and 0.05 at 2096 
/im. 

" Based on surface density diagram in Binney & Merrifield (1998). 
Prom Westerhout (1958). Effective electron density ne// = 10 cm~^ . 
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Fig. 9. Global SED for our Galaxy obtained by summing up the contributions from dust associated with atomic, 
molecular and ionized gas (black solid line), a and T parameters are from Table 4. Also shown are the individual 
emission components: dust in atomic (blue line), molecular (dark green line), diffuse ionized (red line) and compact 
ionized (magenta line). 



